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Abstract 

O" 1 

The Bianchi type III dark energy models with constant deceleration pa- 
rameter is investigated. The equation of state parameter w is found to be 
time dependent and its existing range for this model is consistent with the 
recent observations of SN la data, SN la data (with CMBR anisotropy) 
C^) ' and galaxy clustering statistics. The physical aspects of the dark energy 

models is discussed. 

Keywords : Bianchi type III Universe, Dark Energy, EoS Parameter 






> 



Q ; 1 Introduction 

O 



X 



Arguments have been put forward that we live in a spatially flat low matter den- 
sity universe which is currently undergoing a period of accelerating expansion. 
If the observational evidence upon which these claims are based are rainforced 
; "*] ■ and strengthened by future experiments, the implications for cosmology will be 

incredible. It could then appear that the cosmological fluid is dominated by 
some sort of fantastic energy density, which has negative pressure and has just 
begun to play significant role today. Recent years have witnessed the emergence 
of the idea of an accelerating universe. Therefore, due to some observational 
results pQ — [4] it is now established that universe is accelerating. This signifies 
a remarkable shift in cosmological research from expanding universe to accel- 
erated expanding universe. Now, the problem lies in detecting an exotic type 
of unknown repulsive force, termed as dark energy. The detection of dark en- 
ergy would be a new clue to an old puzzle: the gravitational effect of the zero 
- point energies of particles and fields. The total with other energies, that are 
close to homogeneous and nearly independent of time, acts as dark energy. The 
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paramount characteristic of the dark energy is a constant or slightly changing 
energy density as the universe expands, but we do not know the nature of dark 
energy very well [5] — [T3] . Dark Energy (DE) has been conventionally character- 
ized by the equation of state (EoS) parameter u = p/p which is not necessarily 
constant. The simplest DE candidate is the vacuum energy (lj — —1), which is 
mathematically equivalent to the cosmological constant (A) . The other conven- 
tional alternatives, which can be described by minimally coupled scalar fields, 
are quintessence (lj > —1), phantom energy (w < —1) and quintom (that can 
across from phantom region to quintessence region) as evolved and have time 
dependent EoS parameter. Some other limits obtained from observational re- 
sults coming from SN la data [13] and SN la data collaborated with CMBR 
anisotropy and galaxy clustering statistics [15] are —1.67 < lj < —0.62 and 
— 1.33 < lj < —0.79 respectively. However, it is not at all obligatory to use a 
constant value of lj. Due to lack of observational evidence in making a distinc- 
tion between constant and variable u>, usually the equation of state parameter is 
considered as a constant [THIHI] with phase wise value —1, 0, —1/3 and +1 for 
vacuum fluid, dust fluid, radiation and stiff fluid dominated universe, respec- 
tively. But in general, oj is a function of time or redshift [H] — [20]. For instance, 
quintessence models involving scalar fields give rise to time dependent EoS pa- 
rameter lj [H] — [23] ■ Also some literature is available on models with varying 
fields, such as cosmological model with variable equation of state parameter in 
Kaluza-Klein metric and wormholcs [25, 26 . In recent years various form of 
time dependent ui have been used for variable A models [27] [28] . Recently Ray 
et al [29] , Akarsu and Kilinc [30] and Yadav et al [31] have studied dark energy 
model with variable EoS parameter. 

Spatially homogeneous and anisotropic cosmological models play a signifi- 
cant role in description of the large scale behaviour of universe. Bianchi type 
III cosmological model in presence of dark energy have been studied in general 
relativity by numerous authors. Lorentz [3 2) has presented a model with dust 
and cosmological constant. Chakraborty and Chakraborty have given a bulk 
viscous cosmological model with variable G and A [33]. Singh et al [33] have 
investigated a model with variable G and A in presence of perfect fluid by as- 
suming a conservation law of energy-momentum tensor. Recently, Tiwari [35] 
has studied a model in presence of perfect fluid and time dependent A with con- 
stant deceleration parameter. Bali and Tinkar [36] have investigated a model 
in the presence of bulk viscous borotropic fluid with variable G and A. Un- 
like Robertson- Walkar metric, Bianchi type III can admit a dark energy that 
yields an anisotropic EoS parameter according to their characteristics. The cos- 
mological data from the large scale structure [37] and type la supernova [3j [4] 
observations-do not rule out the possibility of anisotropic dark energy either 
[3J [3J. 

In this paper, we have investigated the anisotropic DE models with variable 
lj. This paper is organized as follows: In section 2, the metric and field equations 
are described. The solution of field equations are presented in section 3 and 
section 4 concludes the findings. 



2 The Metric and Field Equations 

We consider Bianchi type III metric in the form 

ds 2 = -dt 2 + A 2 dx 2 + B 2 e 2ax dy 2 + C 2 dz 2 (1) 

where A, B and C are the function of t only. 

The simplest generalisation of EoS parameter of perfect fluid may be to de- 
termine the EoS parameter separately on each spatial axis by preserving the 
diagonal form of the energy-momentum tensor in a consistent way with the 
considered metric. Thus, the energy momentum tensor of fluid is taken as 

T°=diag[7$,T},I*,li\ (2) 

Then, one may parametrize it as follows, 

Tj = diag [p, -p x , -p y , -p z ] = diag [1, -u x , -U V) -u z ] p 

= diag [1, — u>, — (w + S) , — (u> + 7)] p (3) 

where p is the energy density of fluid,; p x , p v and p z are the pressures and 10 x , uj y 
and lo z are the directional EoS parameters along the x, y and z axes respectively. 
u is the derivation-free EoS parameter of the fluid. We have parametrized 
the deviation from isotropy by setting uj x = uj and then introducing skewness 
parameter 5 and 7 that are deviation from uj along y and z axis respectively. 
The Einstein field equations, in gravitational units (87rG = 1 and c = 1), are 

Rij - \Rgn - ~Tij (4) 

where the symbols have their usual meaning. 

In a comoving co-ordinate system, Einstein's field equation (J4J, for the anisotropic 
Bianchi type III metric ([I}, in case of ((3]), lead to the following system of equa- 
tions 
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Here, the sub indices 4 in A, B, C and elsewhere denote differentiation with 

respect to t. 

Integrating equation ([9]), we obtain 

A = kB (10) 



where k is the positive constant of integration. We substitute the value of 
equation (TTOj) in equation (J7J and subtract the result from equation ([6]), we 
obtain that the skewness parameter on y-axis is null i.e. 

5 = 

Thus system of equations from (0 - (j9|) may be reduce to 

A* A4C4 a 2 

T^ 2 ^c-T 2=p (11) 

-f + ^ + ^ = . up (12) 

2^ + |-^=-(- + 7)p (13) 

Now we have three linearly independent equations (jlip - ()13[) and hve unknown 
parameters (A, B,u,p, 7). Two additional constraint relating these parameters 
are required to obtain explicit solutions of the system. 

(i) The law of variation of Hubble's parameter that yields a constant value of 
deceleration parameter. Such type of relation have already been considered 
by Berman [40] for solving FRW models. Later on many authors (Singh et al 
[4T| [42] and references therein) have studied flat FRW and Bianchi type mod- 
els by using the special law of Hubble parameter that yields constant value of 
deceleration parameter. 

(ii) We assume that the expansion (9) in the model is proportional to the shear 
(a). This condition leads to 

B = C" (14) 

where n is proportionality constant. The motive behind assuming condition (ii) 
is explained with reference to Throne (1967) [43], the observations of the veloc- 
ity red-shift relation for extragalatic sources suggest that Hubble expansion of 
the universe is isotropy today within « 30 percent (Kantowski and Sachs |44j : 
Kristian and Sachs [45). To put more precisely, red-shift studies place the limit 

^<0.3 

on the ratio of shear a to Hubble constant H in the neighbourhood of our galaxy 
today. Collin et al [35] have pointed out that for spatially homgeneous metric, 
the normal congruence to the homogeneous expansion satisfies that the condi- 
tion § is contant. 

The average scale factor of Bianchi type III metric is given by 

R=(ABCe ax Y (15) 



We define, the generalised mean Hubble's parameter H as 

H = i (#1 + H 2 + H 3 ) (16) 

where Hi = 4 1 , H 2 = Mf and H 3 = ^=f are the directional Hubble's parameter 
in the direction of x, y and z respectively. 
Equation (|T5|l . may be reduces to 

Since the line-element ([T]) is completely characterized by Hubble's parameter H. 
Therefore, let us consider that mean Hubble parameter H is related to average 
scale factor R by following relation 

H = kiR~ s (18) 

where k\ > and s > 0, are constant. 
The deceleration parameter is defined as 

RRaa ,.,„n 

From equations (fT7j) and (|18|) . we get 

R A = k x R- s+1 (20) 

r 44 = -fc2( s _ i)ir 2s+1 (21) 

Using equations (p~5]) and (fT9|). equation ([TTll leads to 

g = s — 1 (constant) (22) 

The sign of q indicates whether the model inflates or not. The positive sign 
of q corresponds to standard decelerating model where as the negative sign of 
q indicates inflation. However the current observations of SN la and CMBR 
favour accelerating models i.e. q < 0. 
From equation ([T9l) , we obtain the law of average scale factor R as 
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(23) 



where c\ and c 2 are the constant of integration. 

From equation (|22|) . for s ^ 0, it is clear that the condition for expansion of 
universe is s > i.e. q + 1 > 0. Therefore for expansion model of universe the 
deceleration parameter (q) should be greater than —1. 



3 Solutions of the Field Equations and Discus- 
sion 

3.1 Case(i): when s^O 

Equations ©, (JHH and ([23) lead to 

B = l(Dt + a)r (24) 

Equation ([14]) and fl2j} lead to 

C = ;i(L»t + ci)™ (25) 

From equations (|T0|) and ((24]), we obtain 

A = L(Dt + a)r (26) 

where c is the constant of integration and I — (cq) 2n + 1 , L — kl, l\ — l~ , r = 

(2n+l)s 
3n 

Thus the Hubble's parameter H, scalar of expansion 6 and shear scalar a are 
given by 

H = , h : (27) 

(Dt + ci) 

g = 3 # = 3fel (28) 
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Equations fl28]) and ([29]) lead to 

<r _ £>(n- 1) 
5 \j3nrki 

Using equations (fTTj) , (1231) and ([2^]) . the energy density of the fluid is obtained 

as 

D 2 (n + 2) a 2 

P = 9 T (Ol ) 

nr 2 {Dt + ciY L 2 {Dt + Cl )~ 

Using equations (TT2l . ([25]) . (|26|) and (|3T|) . the equation of state parameter w is 
obtained as 



n 2 r 2 (Dt+ Ci y ,„ 2 < 



_D 2 Q+2) 



^f nr 2 (Di+ci) 2 



L 2 (nt+ci 

Using equations (fT3|) . ([26]). (f3T|) and ([32]) . the skewness parameter 7 (i.e. devi- 
ation from uj along z-axis) is given by 

D 2 [n 2 (2-r)+n(r-l)-l] 
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Figure 1: The plot of EoS parameter (w) vs. time (t) with D = I, n = 0.25, 
s = 0.5 d = .05, a = 0.3 and L = 1.5 



From equation (|32p , it is observed that the equation of state parameter uj is time 
dependent, it can be function of redshift z or scale factor R as well. The red- 
shift dependence of uj can be linear like uj(z) = ujq + uj 1 z with uj 1 — {^) = 
[13 US] or nonlinear as oj(z) = uj + f^f [HI |3D]. The SN la data suggests 
that —1.67 < oj < —0.62 [T3] while the limit imposed on u by a combina- 
tion of SN la data (with CMB anisotropy) and galaxy clustering statistics is 
— 1.33 < uj < —0.79 15! . So, if the present work is compared with experimental 
results mentioned above then, one can conclude that the limit of uj provided by 
equation (|32l) may accommodated with the acceptable range of EoS parameter. 

Also we see that for s — , , -.w" . ;. , the uj vanishes. Thus for this particular 

(n+l)(2n+l) ' f 

value of s, our model represents dusty universe. 



The variation of equation of state parameter (w) with cosmic time (t) is 
clearly shown in Fig. 1, as a representative case with appropriate choice of con- 
stants of integration and other physical parameters using reasonably well known 
situations. From Fig. 1, we conclude that in early stage, the EoS parameter 
uj was positive (i.e, the universe was matter dominated) and at late time it is 
evolving with negative value (i.e. at the present time). The earlier real matter 
later on converted to the dark energy dominated phase of universe. 



In absence of any curvature, matter energy density (f2 m ) and dark energy 



density (Oa) & r e related by the equation 

n m + n A = i (34) 

where J7 m = jjp- and Q& = jjj? 
Thus equation ([34]) . reduce to 

Using equations (j27j) and (J3TJ) . in equation (|35|) . the cosmological constant is 
obtained as 
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(36) 



Measurements of redshift of la supernova SN 1997 ff indicate the present ac- 
celeration of universe. Recently, Carmeli and Kuzmenko [51 have shown that 
the cosmological relativistic theory predicts the value for cosmological constant 
A = 1.934 x 10~ 35 s -2 . This value of A is excellent in agreement with the 
measurements recently obtained by the High-z Supernovae Cosmology Project. 
Recent supernovae la [31 H] observations suggest that the cosmological constant 
is decreasing function of time and they approach to small positive values as 
time increases (i.e. the present epoch). From equation (|36[) . we observe that 
cosmological constant (A) decreases as time increases. 

3.2 Case(ii): when s = 

Equations ©, (dU, flT7J} and (J23J lead to 

B = l e k2t (37) 

Equations dHJ and ([37]) lead to 

C = he^ (38) 

From equations (JTUJ) and (|37p . we obtain 

A = L e k2t (39) 

where /o is constant of integration and Lq = klo, I2 = Iq •, &2 = 2n+\ ■ 

Thus the Hubble's parameter H, scalar of expansion 9 and shear scalar a are 

given by 

H = ( 2 " + 1 >*> (40) 

6n 

a _(2n+l)k 2 (41) 
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Figure 2: The plot of EoS parameter (w) vs. time (t) with a — 0.3, n 
0.25, k 2 = 0.5 and L = 0.10 



Equations flU} and (g^J) lead to 

a (n-1) 



v^n + l) 



(43) 



Using equations (fTTj) . (|3"5)) and (13T)1) . the energy density of the fluid is obtained 
as 

p = ~ — ~^iq (44) 

Using equations (|T2"j). (|55|) . (J3T))) and (|4"4")l . the equation of state parameter a; is 
obtained as 

!-2 

(45) 
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Using equations (fl~3|) . (j39|) , (l44l and (1451) . the skewness parameter 7 (i.e. devi- 
ation from uj along z-axis) is given by 
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The Eos parameter u> is found to be negative and its time varying nature is 
clearly shown in Fig. 2. 



Using equations (|4"0")) and (|44"]) . m equation (j3"5|) . the cosmological constant is 

obtained as 

_ 3(2n + l) 2 K 2 _ (n + 2)k 2 L 2 -na 2 e- 2k * t 

A ~ rw.2 „r2 i 47 ) 
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From equation (|47p . it is clear that cosmological constant A is decreasing func- 
tion of time and approaches to small positive value at late time which is sup- 
ported by results from supernova observations recently obtained by High-z Su- 
pernova Team and Supernova cosmological project [3] 0] • 

4 Concluding Remarks 

In this paper, we have studied anisotropic DE models with variable EoS param- 
eter u, considering two cases, (3.1) and (3.2) for s ^ and s = respectively. 
It is observed that in both cases, EoS parameter uj is variable function of time 
which has been supported by recent observations [UJ [TS] and work of several 
authors (for example, see Refs. [21] — [31]). In case (i), we have shown that, 
in early stage, the equation of state parameter u was positive i.e. the universe 
was matter dominated and at late time, it is evolving with negative values i.e. 
present epoch (Fig.l). Where as in case (ii), for s — 0, from equation (|45]) . we 
have obtained that, at cosmic time t = -^-In-, — " ." ■> , s , lo — — 1 (i.e. cosmo- 

logical constant dominated universe), when t < -^-In -, jfrhzr? ; w > — 1- (i- e - 

quintessence), and when t > ^~l n f~r^Tp ; w < ~ '!• (i- e - super quintessence 
or phantom fluid dominated universe) [52], representing the different phases 
of universe through out the evolving process. Since in both cases, ? = con- 
stant, the models do not approach isotropy at any time. Therefore, we can not 
rule out the passibility of anisotropic nature of DE at least in Bianchi type III 
framework. 
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